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ABSTRACT 


The thesis presents a theoretical investigation carried out 
to study the transient behaviour of an element of a 
thermoelectric cooling device. A mathematical model of a simple 
thermoelectric device has been developed and the resulting 
equations solved under various boundary conditions. The hot and 
cold junction temperatures have been plotted against time for the 
various conditions. In order to carry out a comparative evaluation 
of the influence of the various parameters viz. current and 
thermal load, a wide range of values have been considered for each 
case. 

It is found that Peltier cooling at the cold junction occurs 
most rapidly, followed by heat conduction from adjacent sections 
of the thermoelements. The process of the transfer of the Joulean 
heat generated within the volume of the thermoelements takes the 
longest time, but exercises an overriding influence, especially at 
higher current values. 
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CHAPTER - 1 


INTRODUCTION 

In 1821 Seebeck[13, reported some experiments to the Prussian 
Academy of Sciences, in which he had discovered thermoelectric 
currents arising in a closed circuit made up of different 
conductors kept at different junction temperatures. Through 
extensive experimentation Seebeck was also able to arrange a great 
variety of solid and liquid metals and their alloys in more or 
less the thermoelectric series that is recognized today. 

Twelve years later Peltier[2], a French watchmaker, published 
some results which showed that he had discovered a second 
thermoelectric effect. When a current is passed through a junctior 
between two different conductors there is absorption or generatior 
of heat depending on the direction of the current. This effect is 
superimposed upon, but quite different from, the Joule resistancf 
heating effect usually associated with the passage of an electrif 
current. Peltier, like Seebeck, did not understand the tru' 

significance of his results, but in 1838, Lenz demonstrated tha 
water could be frozen at a bismuth-antimony junction by th 
passage of a current; on reversing the current the ice could b 
melted. 

Thomson[3], later Lord Kelvin, realized that a relatio 
should exist between the Seebeck and Peltier effects and proceede 
to derive this relation from thermodynamical arguments. This l€ 



him to the conclusion that there must be a third thermoelectric 
effect now called the Thomson effect; this is a heating or cooling 
effect in a homogeneous conductor when an electric current passes 
in the direction of a temperature gradient. 

The basic theory of thermoelectricity was first derived 
satisfactorily by A1 tenkirch [4 , 5] , in 1909 and 1911. He showed 
that for both thermoelectric power generation and cooling, 
materials were required with high thermoelectric coefficients, 
high electrical conductivities to minimize Joule heating and low 
thermal conductivities to reduce heat transfer losses. However it 
was quite a different matter knowing the favourable properties and 
obtaining materials embodying them, and so long as metallic 
thermocouples were employed no real progress was made. It is only 
since semiconductor thermocouples have been developed that 
reasonably efficient thermoelectric devices have become possible. 

1.1 Steady State versus Transient Behaviour : 

Steady State Behaviour : The steady state behaviour refers 
to the operation of the thermoelectric cooler under static 
conditions of heat load, current, ambient temperature etc. The 
steady state behaviour of a thermoelectric device has been 
extensively studied and the results can be obtained from any 
standard text on the sub ject [6 , 7 , 8 3 . In particular the effect of 
various thermoelectric parameters on the performance of a 
thermoelectric device is given by the figure of merit, z, which is 
defined as : 


z 


2 O 

a . -t: 



where a is the Seebeck coefficient (in V/K), of the junction and 
is equal to the difference between the two absolute Seebeck 
coefficients of the materials that make up the junction. 
o is the electrical conductivity in 
k is the thermal conductivity in W/m-K 

Under steady state operating conditions the coefficient of 
performance and the maximum temperature drop depend only on the 
value of z, which thus describes the quality of the thermoelectric 
refrigerator element. 

Transient Behaviour : The transient behaviour refers to the 
operation of the device in the initial period, before steady state 
is reached. There is at present very little information available 
about the behaviour of thermoelectric devices under non-steady 
state conditions. This problem is much more difficult than the 
steady state one as it requires the solution of partial 
differential equations that are at best linear in temperature but 
non-linear in current. Further these equations are subject to 
boundary conditions that include temperature-current product type 


non- li near i ties . 




Fig. 2.1. Interdependence of Heat Flux and Electric Current 

The two ends of a metallic bar are in contact with two heal 
reservoirs maintained at different temperatures. If > Tg , the: 
heat is transferred by conduction from reservoir A to reservoir B 
In addition, any mobile charge carriers present, possess a highe: 
average thermal energy at the hot end A as compared to the coo 
end B . Consequently the average velocity component along th 





positive X direction of carriers at A exceeds the average 
velocity component along the negative x direction of the carriers 
at B . This imbalance leads to a net flux of mobile particles 
from A to B until the charges accumulating at B and the resulting 
deficit at A give rise to an internal electric field which 
ultimately stops any further net transfer of carriers. Under 
steady state conditions the concentration of mobile carriers at 
the cold end exceeds that at the hot end, and a difference of 
potential is developed across the bar. During the establishment of 
the steady state, an electric current accompanies the normal heat 
flux. 

Now suppose an external source causes the potential 
difference between the ends of the bar to be altered in such a way 
as to induce a flow of carriers from left to right in the positive 
X direction. Carriers originally at temperature find 
themselves displaced to cooler surroundings, and they discharge 
their excess kinetic energy by collision with the lattice. This 
process assists the normal conduction of heat from the hot to the 
cold junction. If however the potential gradient is so altered as 
to cause a flow of carriers along the negative x direction, then 
particles at temperature Tg accuiriulate at the hotter end, thereby 
altering the heat flux which would be observed in the absence of a 
current . 

2.2 Origin of Thermoelectric Cooling in a Semiconductor [6, 10] : 

A semiconductor may be defined as an electrical conductor 
with a finite forbidden energy gap between the valence and the 
conduction bands. In an insulator the valence band is completely 



full and the conduction band is aipty of electrons. Quantum 
mechanical considerations shew that, in these circunetances, no 
electric current can flow. In an intrinsic semiconductor, charge 
carriers are produced by thermal excitaticn of electrons frem the 
valence band to the conduct im band, across the forbidden energy 
gap. In an extrinsic semiconductor the presence of certain 
inpurities (donors) nay provide electrons in levels which are 
close enough to the conduction band for them to be excited into 
the latter at relatively low tenperatures. Alternatively, other 
impurities (acceptors) may provide vacant levels, close to the 
valence band, into which electrcffis from the latter may pass. The 
vacant spaces in the otherwise full valence band behave like free 
electrons with a positive charge and are usually known as positive 
holes. Ssni conductors are called n - type when the conduction of 
electricity is prinarily due to electrons, and p - type wh^ 
conducticHi is mainly due to positive holes. 

In an electric circuit vhen an electron flows from the 
connecting metal strip into the n - type semiconductor, it has to 
rise into the conduction band of the latter, absorbing an 
equivalent amount of potential siergy in the process. This ^ergy 
absorption at the junction is the origin of the Peltier cooling 
effect, and the Peltier coefficient is, thus, the average value of 
this energy change per unit charge. By convention it is negative 
for an n - type s«ifti conductor and positive for a p - type 


material . 



2.3 Thermoelectric Effects ; 

Three thermoelectric effects govern the behaviour of any 
thermoelectric device in general. These are ; 

(i) The Seebeck Effect 

(ii) The Peltier Effect 

(iii) The Thcmson Effect. 

These effects are a manifestation of the direct reversible 
conversion of electrical energy into thermal energy and 
vice-versa. 


2.3.1 Seebeck Effect : 

If a closed circuit is made of two dissimilar metals, an 
electric current flows in the circuit v^nen the two junctions are 
maintained at different teiperatures . For a aiall tenperature 
differaice betwe^i the two juncticms the open circuit voltage is 
proportional to the tenperature difference t 

AV = a./liT 

vhere AV = open circuit voltage developed 
a = Seebeck Coefficient 

AT = taiperature difference between the two 
junctions. 

For a small taiperature difference AT, this can be writtai as ; 


a 


lim 

AT^ 



dV 

dr 


...( 2 . 1 ) 


which gives us the definition of the Seebeck Coefficient. 


The Seebeck Coefficient depends upon the choice of the two 
materials and cannot be attributed to either material alone. In 



practice, the absolute Seebeck coefficient of a material is 
determined with respect to a material such as lead in which the 
Seebeck coefficient is negligible. The Seebeck coefficient of a 
thermocouple made of two different materials is the difference 
between the absolute Seebeck coefficients of the two materials. 

2.3.2 Peltier Effect : 

Wien a direct current is passed through a junction of two 
dissimilar metals, the junction becomes either hot or cold. The 
rate of liberation or absorption of heat is proporticnal to the 
current and is given by: 

Q = n.I . . .(2.2) 

vAiere Q = rate of heat liberated or absorbed 
n = Peltier Coefficisit 
I = the direct current in the circuit. 

As in the case of the Se^^eck coefficient the Peltier 
coefficient too cannot be ascribed to either material alone but 
rather to the junction as a whole. 

2.3.3 Thomson Effect : 

Unlike the Seebeck and Peltier effects vbich take place caily 
at the junctions, the Thcmscan effect involves the interchange of 
energy with the surroundings throughout the bulk of the conductor 
if an electric current flows through it, while a terperature 
gradient exists across the conductor. The power absorbed or 
liberated per unit length is given by the expression ; 




I = the current through the conductor 
dT 

= tenperature gradient across the ccnductor. 


Again unlike the Seebeck and Peltier coefficients the llicnBon 
coefficient depends solely on the properties of the single 
material . 

The Seebeck, Peltier and Thcmson effects are reversible in 
nature. In additicm to the above three reversible effects 
associated with thermoelectricity, there are two other 
irreversible effects that take place - Jouleean Heating and 
Thermal Conduction from the hot juncticm to the cold junction. All 
the reversible and irreversible effects have to be takei into 
account in thermoelectric cooling analysis. 
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CHAPTER - 3 

LITERATURE SURVEY AND THE PRESENT WORK 

3.1 Literature Survey : 

Though the principles behind thermoelectric cooling have been 
known for a long time, practical devices based on them could only 
be developed recently, after the development of new semiconductor 
materials. Thus thermoelectric cooling is still one of the young 
branches of modern technology. It has assumed an added 
significance in recent years because of the need for elimination 
of CFCs and the consequent need to develops alternative 
refrigeration technologies. 

E. A. Kolenko [11] has described some of the earliest 
engineering designs of thermoelectric cooling devices, developed 
in the sixties at the Institute of Semiconductors of the U.S.S.R. 
Academy of Sciences- However these devices are all based on steady 
state principles only. 

A. F. Ioffe [12], in his pioneering work on thermoelectricity 
has briefly written about the operation of thermoelectric devices 
in the initial stages, before steady state is reached. 

L. S. Stil'bans and N. A. Fedorovich [133 analysed th€ 
operation of refrigerating thermoelectric elements in non-steady 
state conditions. They modelled the cold junction of the 
thermocouple as a lumped mass, and obtained an exponentially 
decreasing temperature profile when the current was passed througl 
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the device. 

A. G. Shcherbina [11] has reported on the operation of 
semiconductor refrigerators in the initial period before steady 
state conditions are reached. In his analysis he has considered a 
single thermoelement and has modelled the heat transfer as a 
volume process. First he has considered the simplest operating 
condition in which the temperature has been taken as constant and 
equal to the temperature of the ambient medium. The second case is 
a modification of the first one, in as much as the hot junctions 
are taken to be in contact with a moving liquid, the temperature 
of which is taken to be constant and equal to the temperature of 
the surrounding medium. In both the cases the final expression for 
the cold junction temperature shows an exponential decrease with 
time. 

A very comprehensive analysis of the transient behaviour of 
thermoelectric devices has been carried out by P. E. Gray [14]. He 
has derived the general governing equations that describe the 
phenomena and has obtained solutions for various operating 
conditions. In his analysis Gray has assumed that each of the 
variables can be written as the sum of a large-signal static or 
quiescent term, which is not a function of time, and a 
small-signal transient or incremental term which is small compared 
to the corresponding quiescent term. The separation leads to 
large-signal static equations and small-signal transient 
equations. The latter are then Fourier analysed to obtain the 
transfer functions which describe the behaviour of the device 
under various conditions of interest. 



3.2 Present Work ; 


In recent years thermoelectric cooling devices have been 
developed for use in nuclear physics, vacuum technology, 
metallurgy, medicine, botany, electronics and many other branches 
of science and technology. Many such devices require the 
temperature to be controlled within a very narrow range. For 
example a thermoelectric spot cooler may be required to control 
the temperature of an electronic component that dissipates heat at 
a variable rate. For the design and development of such devices 
their transient behaviour must first be investigated. 

From the literature survey it is apparent that 
thermoelectricity and its application to refrigeration is still a 
relatively new field. This is especially true of the transient 
behaviour of such devices where the work done mainly involves 
gross approximations, and a rigorous analysis of the phenomena has 
not even been attempted. Hence the need for the present work to 
formulate a mathematical model of an element of a thermoelectric 
cooler. The two thermoelectric effects (Peltier cooling and 
Thomson effect) and the associated effects (Joulean heat 
generation and heat conduction) have been considered in their 
entirety, to derive the generalized governing equations for the 
system. The equations have been solved numerically under various 
boundary conditions of practical interest. Reasonable assumptions 
have been made only where necessary to ensure accurate results. 



CHAPTER - 4 


DEEOVMION CF THE GENERAL GOVEfStlNG EQUATIONS 


4.1 Mathematical Model of the System : 

A thermoelectric cooling element in its sirrplest form is a 
sirrple p-n thenroelectric couple. It consists of two homogeneous 
semiconductors to v^ich electrical connections are made by netal 
strips that are assxjmed to have negligibly small specific heats, 
thermal and electrical resistances . The geometrical arrangement is 
shown scheratical ly in Fig. 4.1 : 


ELB4ENT 1 
T^(x,t) 


ELEMENT 2 

T2(x,t) 


Model of the Thermoelectric Device 





The thermoelernents, which are uniform in cross section are 
assmed to have the same length. The sides of the thermoelements 
are further assumed to be insulated so that the flow of heat is 
along the length of the thermoelements only (one dimensional flow 
of heat ) . 

4.2 Sign ConvKitions Followed : 

(i) The position coordinate, x, is measured downwards frcxn the 
upper junction vhich is ta]<en at x = 0. 

(ii) The heat entering the junctions both at x = 0 and at x = 1 

taken as positive. 

(iii) The current is defined as positive in the positive 
x-direction in element 1 and in the negative x-direction in 
element 2. 

(iv) The terminal voltage is defined such that the rate at viiich 
aiergy leaves the system at the electrical terminals is given by 
V*I watts. 

(v) A positive Thcmscaa coefficient denotes liberation of heat if 

the tenfperature gradient and electric current are in the same 

direction. The net Thcnson coefficient for the couple is given by: 

X = X - X ...(4.1) 

1 2 

(vi) The Peltier coefficient vdiich is defined such that a 

positive coefficient designates heat liberation the current 

flew is positive ( that is from element 2 to elanent 1). Similarly 
the Peltier coefficient is defined such that a positive 

coefficient designates heat absorpticx} for a positive current 


is 


flow. 



4.3 Derivation of the General Governing Equations ; 

The general governing equaticns that ccsrpletely describe the 
behaviour of the cooling elsment shown in Fig. 4.1, consist of two 
differential equations and three energy balance equations : 

(i) Two differential equations that must be satisfied by the 
terperature distributions in the elements. 

(ii) Energy balance for each iiincticn which relates P and P , to 

0 1 

the electric currents and temperatures of the two elerents. 

(iii) Overall energy balance equation. 


4.3.1 Differaitial Equations for Tenperature Distribution : 

Ccsisider a differential strip of elesnent 1 , of length dx 
at a distance x frcxn the origin : 



X 


X + dx 


Fig. 4.2. Differential Strip of Elemeit 1 


The rate at which heat is conducted into the segment at 
X + dx is ; 


P 

1 n 


= 


dr. 


dx 


. . .(4.2) 


x+dx 



ll'ie rate at v^ich heat is conducted out of the segnnent at 


X IS : 


P = k, .a, . 

out 1 1 


ar. 


dx 


...{ 4 . 3 ) 


The rate at vhich Joule heating takes place in the segment is 

Pj = (current) x( resistance of the segment) 
dx 


= I'^.P 


1 a 


... ( 4 . 4 ) 


The rate at v\hich heat is liberated due to the Thcmscan 
effect is : 

dr. 




.dx 


3x 


...( 4 . 5 ) 


The rate at v^iich stored energy is increasing is : 

ar. 


P = c .a . 

s 1 1 


.dx 


...( 4 . 6 ) 


dx 


Conservaticxi of energy demands that : 

Net energy conducted into the segment + Joule heating within the 
segment + Heat liberated due to the Thcrscxi effect = Rate at vhich 
the stored energy increases. 


ie. 

p 

1 Ti 

- p 

out 



= P 

or 

P. 

1 n 

~ p . 

out 

= p 

8 

- Pj 

- P, 



le. 




ar. 


dx 


“ k .a . 
1 1 


x+dx 




9x 


c .a . 
1 1 


ar 


9x 


.dx 


ar. 


t. 2 ^ dx ^ 

- I .p . - I.'C . 

la ^ ^ 

1 ax 


.dx 


le. 


a 

ax 


ar^ >1 ar 

k,.a,.— J.dx = c,.a,.— - l".p 

^ ax '' dy. 


dx 

a. 


I.r^. 


3T. 


dx 


.dx 


le. 


dx 


k,.a,. 


ax 


Pi 2 

+ t .1. + — .1 - c,.a . — = 0 

1 ^ a 11^ 

ax 1 dx 


...(4.8) 


is the differential equatiOTi that governs the teirperature 
distribution in elanent 1. The derivation of the differential 
equation for elanait 2 is identical except that the current I is 
in the opposite direction and so must be r^laced by -I. Haice the 
equation for the tsrperature distribution in element 2 is : 



.a 


2‘ 



t3.I- 


ar. 


dy. 




= 0 


...(4.9) 


4.3.2 Jiincticai Energy Balance : 

This is obtained by recognizing that the power vA^ich enters 
either the upper or lower junctions is the sum of the powers that 
are rsnoved from the corresponding junctions by the Peltier effect 
and by thenral conduction through the the tvK> elements. Thus, 



OT {x,t) ar (x,t) 

P (t) = -n .I(t) - k .a . - k,.a,. — ^ 

O O 11- 22^ 

OX - OX „ 

x=0 x=0 


...(4.10) 


ar (x,t) aT,{x,t> 

P (t) = n .I(t) + k -a . +k,.a,. — = 

i £ i i ^ 

OX , dx , 

X=i X=1 


Now the first Kelvin relation gives us : 


. . . (4.11) 




. . (4.12) 


a .T = n 
o o o 


. . (4.13) 


and the second Kelvin relation gives us 


o _ _T_ 

dr T 

o o 


...( 4 . 14 ) 


1 _ 


..( 4 . 15 ) 


[ Please see Appendix - A for the derivation of the Kelvin 
relations.] 

Using equations (4.10) and (4.11) we get : 


- h-*.- 


aT^(x,t) 


k, .a, . 
2 2 




...( 4 . 16 } 
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and 


Pj (t) 


.T|t) .I{t) + 


ax^ (x, t ) 
dx 


+ k . a „ . 


ax^ (x, t ) 


X=1 


ax 


X=1 

..{4.17: 


which represent the energy balance at the two junctions 


4.3.3 Overall Energy Balance Equation : 

Xhe difference between the net thermal power input and the 
electrical power output must equal the rate at which the stored 
energy in the system is increasing : 


ie. P(t)+P^{t)-V(t).I{t) 
1 o 


dW(t) 

dt 


...(4.18) 


where W(t) denotes the total energy stored in the system, 

which is simply the energy accounted for by the specific heats 

of the elements : 

1 1 

W(t) = a . f c .X (x,t).dx + a, . f c^ . X, ( x , t ) . dx ...(4.19) 

ajii ^ j z z 

0 0 

Substituting for W(t),^P (t) and ^P (t) into equation 4.18 
and rearranging the terms we get ; 


V 


f a .T 

- a .X 1 


[(t - T ) .dx 

1 i i 

0 

0 


J 1 2 


i i 

* ^ rs 


. . .(4.20) 



Since the variations in heat load and consequently current 
are generally small all the parameters except resistivity may be 
assumed to be independent of tarperature. 

The reasons for this are : 

(i) Resistivity shows the greatest variation amongst all the 
parameters. 

(ii) Resistivity enters the equations as the coefficient of the 
Joule heating term in the differential equations and this term has 
a quadratic current dependence. 

Further the resistivity can be assumed to have a linear 
variaticn with temperature [7,8] : 


ie. p = p + n .T 

^10 'i 1 

P = P + ri .T 
^2 ^20 *2 2 


Subst i tut i ng these express i ons 


...{4.21) 
. . .(4.22) 

for resistivities into 


equations 4.8 and 4.9 we get s 


a"T, ar ar n 

k .a . + r .1. - c .a . + .1 .T, 

11^2 1 11 _ a^ 1 

dx ax at 1 


+ .1^ = 0 ...(4.23) 

^1 


2 2 


a^ 

2 

T 

T" T 

^2 

r* A 

ax^ 

dx 

at 


% 2 

+ — .1 .T, 
a 2 

2 



= 0 ...{4.24) 



Also equation 4.20 siirplifies to : 


V = f a .T - a -T ] - (T - t ) . (T - T ) 

111 00 / 1 I o 


1 1 

" T -Ko-* 

'■ 1 rs 2 ^ 


1 1 

- + '■'rh*') + i^<p 2 o-i + '■' 2 -Kdx)_ 


V = a 


<-T-a.T l-(t -r).(T-T) 

i i O O J 1 21 O 




-•I 

i. J 1 


T., dx 
2 


...( 4 . 25 ) 

Thus the ccnplete system equations that describe the generic 
thermoelectric device are : 

a"T sr ar r, 

k .a . -f z . 1 . - c ,a . 4 - — - .1 .T 

^ ^ ax" ^ ax ^ ' at ^ 

+ .l" = 0 ...( 4 . 23 ) 


k^.a, . 



T 

T' T 

^2 

ax" 

L • JL * 

" ax 

d _ * 

^ at 


^2 " 


+ .l" = 0 

^2 


...( 4 . 24 ) 



’‘o"> - 


ar,(x,t) 

dx 


- k, .a, . 
2 2 


x=0 


dri^{x,t] 

dx 


x=0 


(4.16) 


Pj(t) = 


a^.T^t) .I(t) + k^.a^. 


aT^{x,t) 

dx 


+ k, .a, . 
2 2 


x=l 


aT2(x,t) 

dx 


x=l 


. . (4.17) 


and lastly. 


V = 

If -T, 

- a .T 1 

- {x - xj 

i.(T - T ) 


[1 1 1 

o o J 

1 2 

1 o _ 


, ( ^10 , ^20 1 .tF r r„ j . ^^2 r j ' 

- - H'57'J -sj-J ^2 j 


(4.25) 


4.4 Alternative Expressicai For The Heat Rejected at the Hot 
Jxmction : 

The rate at which heat leaves the juncticn at x = 1 is also 
given by : 


Rate at which 
heat leaves 
at X = 1 


Rate at which 
heat enters + 
at X = 0 


Rate at vhich 
Joule heat is 
generated 


Rate at vA\ich 
- heat is 
stored 


and following the proper sign convention we have : 


- P (t) = P^(t) - V.I 
1 o 


aw(t) 

at 



ie. - P (t) = - a .T .1 - k -a . 

1 o o 11 


dx 


- ^2-^; 


x=0 


ar. 


5x 


x=0 


] - I's - 

j 1 2 




1 


ar ar 

or P, (t) = - a .T .1 + k .a . + k .a . — ^ 

1 II ^ 1 rs. rt 2 2 

ax x=0 dx 


x=0 


- I 


(r, - t^).(T^ - T^) 




^ ^ r\ ys ^ 


...(4.26) 


4.5 Boundary and Initial Ccaiditicais : 

The teiperature profile in the two thermoelernents is 
determined fay equations (4.16), (4.17), (4.23) and (4.24), subject 
to an initial condition and four boundary conditions. 

The initial condition is given by the fact that the entire 
systan is at the ambient taiperature initially \*ben the current 
flow is nil. Thus : 

T = T, = T for all X ^^ben t 0 

1 2 A 

\faexe T is the ambient temperature. 

A 


dx 



Two boundary conditions follow frcm the physical requirement 

that the terrperatures at each of the two junctions imist be the 

same for both the thermoelements. therefore we have ; 

T = T at X = 0 and at x = 1 
1 z 

The other two boundary conditions are specified by the 
thermal load and current flow through the couple. The different 
situations considered in this work are as outlined below : 

(i) In the first case the effect of a change in the control 

current, on the junction terperatures is analysed. The power input 
at the cold junction, is taken to be constant v^hile the 

power rejected at the hot junction, P^{t), is evaluated frcm 
Eq.4.26. Two cases are have been considered under the inposition 
of these thermal loads, namely : 

(a) there is a step change in the current, and 

(b) the current varies sinusoidally. 

(ii) In the second situation considered the curr^t flow is 
kept ccanstant while the thermal load is varied, to study its 
effect on the junction tsrperatures . Again two cases have been 
considered : 

(a) there is a step change in the power input, P^Ct), and 

(b) the power input varies sinusoidally. 

The power rejected at the hot junction, P^{t), is evaluated frcm 
equation 4.26 in each case. 

The discretisation of the differential equations and the 
method for their solution is given in App^^ix-B. The discussion 
of the results follows in the next chapter. 
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CHAPTER - 5 

RESULTS, DISCUSSIONS AND CONCLUSIONS 

5.1 Thermoelectric Parameters of the Device ; 

The thermoelectric device that has been analysed in this work 
is a p-n couple of bismuth telluride. The p-type element is the 
pure crystal , whi le the n-type element is bismuth telluride doped 

“ 4 

with 5 X 10 percent of iodine by weight. 

The pertinent thermoelectric properties of the elements as 

obtained from the various data plots of bismuth telluride 

[14,15,16,17] are given below : 

k =2.4 N/mK 
1 

k, = 2.2 W/mK 
2 

c = 1.23x10^ J/m^K 

1 

c^ = 1.23x10^ J/m^K 

In the expression for resistivities : 

P, = ... (4.21) 

and 

p = p, + n .T, ...(4.22) 

^2 20 '2 2 

we have ; 

- 5 

p = 1.3x10 ohm-m 

^10 

p = 1.3x10'^ ohm-m 

2 0 

r/^ = 0.814x10"^ ohm/K 

n = 0.814x10"^ ohm/K 



The geometrical parameters of the device are taken to be : 

2 

= 0.836 cm 

2 

a = 0.836 cm 
2 

1=2.0 cm 

5.2 Results and Discussions : 

CASE - i(a) : Step change in the current, with the power input at 
the cold junction kept constant at 0.5 W. 

Thirteen values of the step current are considered from 0.5 A 
to 10 A, to illustrate how the cold junction and hot junction 
temperatures vary with increasing current. The incremental cold 
and hot junction temperatures, ie. (T - T 1 and (T, - T.) , where 

O A 1 A 

T is the ambient (initial) temperature, have been plotted against 
time for the various current values in Figs. 5.1 through 5.11. 

Pigs. 5.1 through 5.4 are for small step current values of 
0.5 A, 0.75 A, 1.0 A and 1.5 A respectively. Initially the hot 
junction temperature is seen to increase while the cold junction 
cools down rapidly. The hot junction temperature reaches a peak 
and then decreases to a final steady state value. 

Two features of the plots are worth noting at this stage. 
Firstly, the initial temperature drop in case of the cold junction 
and temperature rise in case of the hot junction, respectively, 
keeps increasing as the current increases. Secondly the final 
steady state value of the incremental cold junction temperature 
decreases while that of the hot junction increases as the current 


is increased. 



The general trends depicted by the plots show how the device 
responds transiently to the three effects that occur when the 
current is switched on, namely the thermoelectric effects (Peltier 
and Thomson cooling), Joulean heat generation and heat conduction. 
The first observation is explained by the fact that the Peltier 
cooling effect, which is proportional to the first power of the 
current, is a surface phenomenon and has a very low inertia, 
whereas the process of transfer of Joulean heat, generated within 
the volume of the thermoelement, to the cold junction, is 
characterized by a much higher inertia. The rate of absorption of 
Thomson heat, and its effect on the cold junction temperature, is 
same as that of the Joulean heat generation. Thus it serves only 
to modify the Joulean heat, and since its effect is very small it 
will not be considered explicitly henceforth. 

The cold junction therefore cools down rapidly at first as a 
result of Peltier cooling, but the Joulean heat generation and 
heat conduction from the hotter neighbouring areas, soon arrest 
any further cooling of the cold junction till it reaches a final 
steady state value. The hot junction temperature shoots up 
initially due to Peltier heating but heat conduction to the 
neighbouring layers soon causes the temperature to drop. The 
Joulean heat generation further modifies the hot junction 
temperature which ultimately approaches a final steady state 
value . 

The steady state temperatures for step currents from 0.5 A to 
1.5 A, are given in Table 5.1. The increasing difference, between 
the hot and the cold junction temperatures, shows how the Peltier 
cooling and heating increases as the current increases. However, 



this is not the complete picture as the Joulean heat generation 
and heat conduction also increase with increased current values. 
Their effects will be evident when we consider higher current 
flows through the device. 

Figs. 5.5 to 5.13 show the behaviour of the thermoelectric 
junction for step currents of 2.0 to 10.0 A. The cold junction 
temperature now shows a very interesting trend. It decreases 
sharply at first, as was seen earlier, but this rapid decrease in 
temperature is soon checked, so that it reaches a minimum value 
and then it increases till it reaches a final steady state value 
much above the minimum. The minimum temperature reached decreases, 
while the final steady state value increases as the current 
increases . 

The hot junction temperatures for step currents of 2.0 A to 
10.0 A (Figs. 5. 5 to 5.13)show a similar trend to Figs. 5.1 to 5.4. 
The steady state temperature increases rapidly as the current 
increases . 

These trends show very clearly the transient action of the 
different effects on the thermoelectric device. The Peltier effect 
occurs most rapidly (almost instantaneously) and so the cold 
junction temperature decreases rapidly as soon as the current is 
switched on. However, heat conduction to the adjacent regions 
arrests the rapid decrease in temperature . Joulean heat 
generation and its transfer to the junction occurs last of all, 
but at higher current values, this is clearly the overriding 
factor. Thus the temperature of the cold junction soon starts 
increasing and for currents of 4.0 A and more the steady state 
temperature is actually higher than the initial temperature which 



means that the device would be quite ineffective as a 
thermoelectric cooler at steady state. 

The hot junction temperatures too illustrate the transient 
action of the three effects clearly. Thus the hot junction 
temperature shoots up initially due to the Peltier heating. Heat 
conduction to the adjacent layers soon follows and is shown most 
clearly in Figs. 5.1 to 5.4 where there is a small but clearly 
defined temperature drop before Joulean heating causes the 
temperature to increase again till a final steady state value is 
reached. The final steady state temperature increases rapidly as 
the current increases since the Joulean effect depends on the 
square of the current. Also the steady state is attained after a 
longer time interval, when the current increases, which further 
shows that Joule heat generation and its transfer have a much 
higher inertia. 

The minimum cold junction temperature and the final steady 
state temperatures of both the hot and the cold junctions, for 
currents of 2.0 A to 10.0 A are also given in Table 5.1. The 
minimum cold junction temperature decreases almost linearly as the 
step current increases, since it represents mainly Peltier 
cooling. However the final steady state values at both the hot and 
the cold junctions increase rapidly as the current increases, 
since it depends mainly on Joulean heating, which in turn has a 
square current dependence. 

The cold junction temperatures for the various step currents 
have been plotted in Fig. 5. 14. The plots clearly illustrate the 
various trends discussed above. Similarly the hot junction 
temperatures have been plotted in Fig. 5. 15 to enable a comparison. 



of the response of the hot junction, to increasing step current 
values . 

In order to illustrate the transient behaviour of the element 
for increasing current values, the incremental temperatures at 
various instants of time have been plotted against current in 
Figs. 5.16 and 5.17. Figs. 5.16(a) to (1) are for the incremental 
cold junction temperatures. Pigs. 5.16(a) to (e) show that 
initially the temperature decreases with an increase in the step 
current. As time increases the temperature at first decreases with 
current till it reaches a minimum value at a particular current. 
For currents larger than this optimum value the cold junction 
temperature shows a rapidly increasing trend as depicted in Figs. 
5.16(f) to (1) . 

The incremental hot junction temperatures have been shown 
plotted against current, at various instants of time, in Figs. 
5.17(a) to (1). The hot junction temperature always increases with 
an increase in the step current. Initially the increase in 
temperature with current is less pronounced as it depends on the 
Peltier effect alone. As time increases the additional heating due 
to the Joulean effect causes the hot junction temperature to 
increase rapidly with the increasing current. 
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Incremental Junction Temperatures for a Step Current 
0.5 A and Constant Power Input of 0.5 watts. 
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Incrannental Jvinction Tannperatures for a Step Current of 
0.75 A and Constant Power Input of 0.5 watts. 




Fig. 5.3 Incremental Junction Tenperatures for a Step Curr^t of 
1.0 A and Ccmstant Power Input of 0.5 v^ratts. 
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Incremental Junction Tannperatures for a Step Current of 
6.0 A and Constant Power Input of 0.5 watts. 
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Fig. 5.10 Incremental Junction Tetiperatures for a Step Current of 
7.0 A and Constant Power Input of 0.5 watts. 
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Fig. 5.11 Incremental Junction Tenperatures for a Step Current of 
8.0 A and Constant Power Input of 0.5 watts. 




Fig. 5.12 Incremental Junction Tenperatures for a Step Current of 
9.0 A and Constant Power Input of 0.5 vyjatts. 



Fig. 5.13 Incrainantal Junction Tamperatures for a Step Current of 
10.0 A and Constant Power Input of 0.5 watts. 
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Fig. 5.14 Cortparison of the Incranantal Cold Junction 
Tenperatures for Increasing Step Current Values. 
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orison of the Incremental Hot Junction 
iratures for Increasing Step Current Values. 
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Fig. 5.17 (a) Incremental Hot Junction Temperatures for 

Increasing Step Current Values , at Various 
Instants of Time . 
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Fig. 5.17 (b) Incrar^ntal Hot Junction Teiperatures 

Increasing Step Current Values , at V 
Instants of Tim© . 
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Fig. 5.17 (c) Incremental Hot Junction Tenperatures 

Increasing Step Current Values , at Vi 
Instants of Time . 








(lASE - i(b) : Sinusoidal change in the currt:.nf 


with the power 


input at the cold junction kept constant at 0.5 w . 

The current is taken to be of the form : 

T = I + I„.sin(«.t) 

0 1 0 

where 1^^ is the fixed component of the current while is the 

amplitude of the variable part. Two values of the current are 
considered, namely : 

I^= 0.5A, and 

I = I = 5.0 A . 

01 0 

Three values of the frequency, w, are taken for each case 
0.1 Hz., 1.0 Hz. and 10.0 Hz. The incremental junction 
temperatures have been plotted against time in Figs. 5. 18 through 
5.23. The current is also plotted above the temperatures in each 
case . 


The figures illustrate the phase behaviour of the device 
under a sinusoidal current input. The hot and the cold junction 
temperature are always in phase with one another and have the same 
frequency as the current. However they have a definite phase 
difference with the current, which increases as the frequency of 
the input current increases. For smaller current values the phase 
difference increases at a slower rate with increasing frequencies. 
The phase difference for the two current values and three 
frequencies are given in Table 5.2. 



The phase difference occurs because the temperature of the 


thermocouples is not a 

1 inear 

function 

of 

the 

current , 

since 

Joulean heat generation 

is a 

function 

of 

the 

square of 

the 

current. This leads to a 

phase 

difference 

between 

the current 

and 


the temperature. The linear components (Peltier effect and heat 
conduction) have a proportionately larger effect at low current 
values, as compared to the Joulean heat and so the phase 
differences too increase at a slower rate for smaller currents. 
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Fig. 5.18 Increnental Junction Teirperatures for a Sinusoidal 
Current of Anplitude 0,5 A eind Frequency o i Hz. with a 
Constant Power Input of 0.5 watts. 



Fig. 5.19 Increinantal Junction Tenperatures for a Sinusoidal 
Current of Anplitude 0.5 A and Frequency l o Hz. with a 
Constant Power Input of 0.5 watts. 




Fig. 5.20 Incremental Junction Tenperatures for a Sinusoidal 
Current of Anplitude 0.5 A and Frequency 10 Hz. with a 
Constant Power Input of 0.5 watts . 
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Fig. 5.21 Incranental Junction Tenperatures for a Sinusoidal 
Current of Anplitude 5.0 A and Frequency 0.1 Hz. with a 
Constant Power Input of 0.5 watts. 




(duuB) T-uajjno 


< — 


(>l) -clujai iB-iuaujeJoui 


Fig. 5.22 Incremental Junction Terrperatures for a Sinusoidal 
Current of Airplitude 5.0 A and Frequency 1.0 Hz. with a 
Constant Power Input of 0.5 watts. 
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Fig. 5.23 Incremental Juncticai Tenperatvires for a Sinusoidal 
Current of Airplitude 5.0 A and Frequency 10 Hz. with a 
Constant Power Input of 0.5 watts. 


TABLE - 5.2 


PHA SE DIFFERENCE BETWEEN CURRENT AND TEMPERATURE 

FOR VARIOUS FREQUENCIES 








rASE - ii(a) : Step change in the power input at the cold 

-junction with the current kept constant at 1.0 A . 

A step power input is given to the system operating initially 
under steady state conditions of current flow of 1.0 A and initial 
power input of 0.5 W. Five values of the power input are 
considered from 0,1 W to 0.5 W. The incremental cold junction and 
hot junction temperatures are shown plotted against time in 
Figs. 5. 24 through 5.28. 

Both the cold and the hot junction temperatures increase as a 
result of the additional heat load. The cold junction temperature 
however increases more rapidly than the hot junction temperature. 
Both the junctions finally attain new steady state values. 

The cold junction heats up almost immediately as a result of 
the additional power input. As the cold junction temperature 
increases the conduction of heat from the hot junction to the cold 
junction decreases and so the hot junction temperature also 
increases. Due to a finite heat capacitance (ie, ability to store 
heat) , of the thermocouple elements, the hot junction temperature 
increases at a slower rate as compared to the cold junction 
temperature. 

As the step power input increases the final steady state 
temperatures of both the cold and hot junctions increase, but the 
temperature difference between the two ends decreases. Thus more 
heat is now being pumped across the thermoelectric couple but 
across a smaller temperature difference. 

The incremental cold and hot junction temperatures for the 
various step power inputs have been plotted in Figs. 5.29 and 5.30 



respectively. Fig. 5.29 shows how the incremental cold junction 
temperature increases as the step power input is increased. 
Similarly Fig. 5. 30 depicts the response of the hot junction to 
increasing values of the step power input. The plots clearly 
illustrate the various trends discussed above. 

The transient response of the thermoelement for increasing 
power values, at various instants of time, has been depicted in 
Figs. 5.31 and 5.32. Figs. 5.31(a) to (1) show that the 
incremental cold junction temperature increases linearly with the 
thermal load at the cold junction. Figs, 5.32(a) to (1) show a 
similar trend for the hot junction temperature. 
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Fig. 5.31 (a) Incremental Cold Jvinction T 

Increasing Step Power Inputs 
Instants of Time . 
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Fig. 5.31 (b) Incaremental Cold Jwiction T 

Increasing Step Power Inputs 
Instants of Time . 
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Incremental Cold Junction T 
Increasing Step Power Inputs 
Instants of Time . 
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Fig. 5.32 (a) Incremental Hot Junction T 

Increasing Step Power Inputs 
Instants of Time . 
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Fig. 5.32 (b) Incremental Hot Jimction Temperatures 

Increasing Step Power Inputs , at Va 
Instants of Time . 
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Fig. 5.32 (c) Incremental Hot Junction Torperatures for 

Increasing Step Power Inputs , at Various 
Instants of Time . 



r change in the power input at t.hp cold 

junct ion witl) the current kept constai^t at 1 . 0 A. 


Tlie syst«T» re^r^nse is considered for a sinusoidal power 
ini)Ut given to a thermocouple initially operating at a constant 
current of 1.0 A and a constant power input of 0.5 W. The 
sinusoidal pc:iwer input is of the form : 


p = p . sin(o. t ) 

wtiere is the arrplitude of the sinusoidal power input and is 
taken to erpial to 0.1 W. Three values of the freqijency, o, are 
considered - viz. 0.1 H 2 . , 1.0 Hz. and 10.0 Hz. The incremental 
cold and hot junction temperatures are shown plotted against time 
in Figs. 5.29 througli 5.31. The power input is also plotted above 
the traiperatures in each case. 

Tlie plots depict the phase behaviour of the device under a 
sinusoidal jrower input. The hot and cold junction teriperatures are 
alwjiys in phase with one another and with the power input. The 
temperatures also have the same frequency v^iich is equal to the 
frerpiency of the sinusoidal power input. The tarperatures are now 
in phase with the sinusoidal power input because the temperature 
is a line^tr function of the power input. 




Fig. 5.33 Incremental Jvinction Tenperattires for a Sinusoidal 
Power Input of Anplitude 0.1 W and Frequency 0.1 Hz. 
with a Constant Current of 1.0 A. 



Fig. 5.34 Incremantal Juncticxi Tenperatures for a Sinusoidal 
Power Input of Anplitude 0.1 W and Frequency 1.0 Hz. 
with a Constant Current of 1.0 A. 



Fig. 5.35 Incremantal Junction Tenperatures for a Sinusoidal 
Power Input of Anplitude 0.1 W and Frequency 10 Hz, 
with a Constant Current of 1.0 A. 


5.3 Conclusions : 

The transient behaviour of a thermoelectric cooling device 
has been thoroughly investigated in the present work. The 
transient response of the device to a change in current is very 
important from the point of view of the controllability of the 
device, as current is the principal controlling parameter in the 
operation of the device. 

The following conclusions can be drawn from the analysis : 

(i) There is an optimum current at which the steady state 
temperature drop of the cold junction is maximum. For currents 
larger than the optimum value the steady state temperature of the 
cold junction increases, instead of decreasing, as the current is 
i ncreased . 

(ii) For currents larger than the optimum value there is an 
initial temperature drop at the cold junction which increases with 
current. This suggests that pulsed current inputs much larger than 
the optimum value could be used to produce temperature drops far 
greater than those achieved with steady currents. The current 
pulses must be timed so as to produce the additional temperature 
drop as soon as the minimum temperature due to the previous pulse 
i s produced . 

(iii) An increase in the thermal load causes the cold junction 
temperature to increase linearly. For currents less than the 
optimum value, an increase in the control current causes a drop in 
the cold junction temperature. But if the current is greater than 
the optimum value then as steady state is reached the cold 
junction temperature increases due to Joulean heat generation. 
Thus the only solution in this case would be to have additional 



cooling elements in parallel, to take up the increased thermal 
1 oad . 

(iv) The cold junction temperature shows a sinusoidal response 
when the thermcil load at the cold junction varies sinusoidally. 
Again for a sinusoidal current the cold junction temperature also 
varies sinusoidally, but there is a phase shift between the 
current and the temperature in this case. Therefore for a 
sinusoidal thermal load the cold junction temperature can be 
controlled by a sinusoidal current if the phase difference is 
taken into account. It must be noted here that in practical 
situations the thermal load will rarely be perfectly sinusoidal. 
However any periodic load can always be written as a Fourier 
series which would amount to the superposition of a number of 
sinusoidal loads. The control current would accordingly have to be 
the sum of a similar Fourier series with the phase difference 
added to each of the terms of the series. 

5.4 Suggestions for Future Work : 

The present study is a step towards the design and 
development of practical thermoelectric cooling devices in which 
temperature controllability is an important factor. In this work 
the temperature response characteristics of only a single element 
have been considered. However thermoelectric devices rarely 
consist of a single stage. Usually two or three stages are 
cascaded to achieve the desired temperature drop. The dynamic 
behaviour of such a system needs to be studied along simil 
lines, as the control aspects of this type of a device are 
interest in the design of practical thermoelectric cooling 



devices . 


Today a large number of semiconductor materials are 
available. Work needs to be done to investigate the effects of 
their property values on the behaviour of thermoelectric devices 
and to find the optimal operating conditions for each material. 

Work is also required to investigate as to how a 
semiconductor element responds to a pulsed current and to find the 
time period of the pulses that would give the largest drop in 
temperature at the cold junction. The control aspects of the 
element under a sinusoidal current also need to be studied 
further, to determine a relationship for the phase difference 
between the cold junction temperature and the current. These 
results would be extremely important for the design of 
thermoelectric cooling devices that have to operate under a 
periodic thermal load. 

Lastly the results established here and in any future work 
need to be verified through extensive experimentation. 
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APPENDIX - A 


DERIVATION OF THE KELVIN RELATIONS 


The Kelvin relations may be derived by applying the laws of 
thermodynamics to the simple circuit shown in the figure below: 



Fig. A. 1 . Simple Thermoelectric Circuit 


A eurrent I passes round the circuit which consists of 
conductors 1 and 2 with .junctions at temperatures T^ and T^ . From 
the principle of conservation of energy the heat generated must be 
equal to the consumption of electrical energy. If the current is 
small enough. .Joule heating may be neglected. Thus, 


,J, ^ * .[*( ^ 

J 




I 

a . I . dT 

J 12 


(A. 1) 


T. 

u 


Differentiating, we get : 
dll 


t 


(A. 2) 


(X 


1 2 


dT 


.•4- 


1 
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In order to obtain a second relation between the coefficients 
it is necessary to apply the second law of thermodynamics to the 
above thermoelectric circuit. In order that this step should be 
valid it is essential that the process should be reversible. 
Though the thermoelectric effects are reversible, but they are 
accompanied by the irreversible processes of Joule heating and 
heat conduction. Thus, the application of reversible 
thermodynamics is not strictly justified. However the use of 
Onsager's reciprocal relations, which are based on irreversible 
thermodynamics, leads to the same concl us i ons [ 1 8] , so we shall 
assume that there is no change of entropy round the circuit of 

Fig. A.l. Then, 


(^) 


- 1 . dT = 0 


(A. 3) 


On differentiating it we get ; 


dn 

12 


+ = 0 


(A. 4) 


Combining equations A. 2 and A. 4 we have : 


12 1 


. . . (A. 5) 


which is the first of the Kelvin relations. It relates the 
Seebeck and the Peltier coefficients. When applied to the two 
junctions of the generic thermoelectric device depicted in Fig. 
4.1, it yields the following equations : 
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and 


u . T - 
0 u 


n 


. . (4.12) 


. . (4.13) 


The second Kelvin relation connects the Seebeck and Thomson 
coefficients. Its derivation is based on the principal of 
conservation of energy. In section 4.3.3 it was shown that the 

overall energy balance for the generic thermoelectric device 

Fig. 4.1 results in the following equation : 



f-i 1 

0 

H 

0 

1 

- J'S - J 

T. 


- I - f-i" 

^ 1 n 0 


0 

The first Kelvih relation, equation A. 5 has been used in the 
derivation of equation 4.20. If the electric current 
have : 

...(A.6) 

Vo= ».-T-»o'"o 

Tq 

Where V„ is the open-oiroult terminal voltage. If equation 
.A. 8 IS differentiated with respect to the result is : 

. . . (A. 7) 


dV . d« 

^ ^ 'Y , 4- oc 

dT, ‘ dT^ 


- ^ 2 ' 


since 


dV 

dT, 


0 


= « ^ by 


definition of the Seebeck coefficient, 


equation A. 7 becomes : 
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d« 


dT 


1 


1 1 
and putting 

couple we get : 


S- "2 


. . . (A. 8) 

V , the net Thomson coefficient for the 


da. 


dT 


X 

T. 


. . . (A. 9) 

1 " 1 

which is the Second Kelvin relation for the junction at 
X = 1. If equation A. 6 is differentiated with respect to T^ and 
if the definition of is used, namely = -CdV^/ dT^ ) , the 

result is : 


da 


dT, 


V 

T. 


. . (A. 10) 


1 ‘ 1 

which is the second Kelvin relation for the junction at x — 0. 

Although equations A. 9 and A. 10 have been obtained by setting 
j = 0 , they are valid when I 5 ^ 0 , since they are relations among 
the parameters of the thermoelectric device, and those parameters 
do not depend on the current. 
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from the above we have 




. . . (B. 1 ) 



rearranging it we have : 


ai. 

k. 

2 

1 . I 

2 

- 

'■'2' 

,1^ 

. T - 

^20-l' 

at 

^2 

dx^ 

4- “ . 

c _ . a . 

Z 2 

ax 

°2 

2 2 
. a . 

2 

2 

^2-^2 
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and on discretising it we get : 


„n + l_ 

^ 2 , i 2 . i 


k, t" . - 2.T" .+ t" ^ 1- .1 - t" 

2 2,1 + 1 2 > i 2 , 1-1 2 2,1 + 1 2 , 1-1 

c, Ax^ o^-SL 2. Ax 

2 ^ 


'h ■ I 

r - T^.i 

‘=2-^2 


^20 ■ ^ 


c . a 
2 2 


= 0 


where , 

denotes the temperature at the i point and the n time 
step in element 2, and so on, 
the other symbols have their usual meanings. 

Rearranging the above we get : 


ie . T, 


^''2 ^2'^ At ..pH ^ 

, " • ^2 , i + 1 

Ax 2.a^ 


2.k n I ] 

1- +— y.At + 1 .t" . 

* 2 2 I 2 , J 


lilil 4L - t" + - At 

2.aJ- ^2^^^ • c,.a^ 


. . (B.2) 


The expression for the power input at x - 0 (cold junction), 


LS derived in section 4.3.2 is : 


P (t) = - « .T(t).I(t) - k.a. 

o O O 11 


aT^(x, t) 


- k_ . a . 
2 2 


ax^Cx, t) 


...(4.16) 


'ain discretising it we have : 


p. - - %-Trr-^ - 


pH + ^ _ t" ^ 

S 2 11 


k ^. a ^. 


n+l _ pn+1 
^22 21 


where denotes the temperature of the first element at the 

location specified by i = 2, and so on for the other terms. 



Now the first boundary condition stipulates that the 
temperatures of the two elements are always equal at x - 0. 

Therefore : 

T = T for al 1 times . 

11 2 1 

Hence the above expression can be rearranged as : 



. I + 
u 




n+ 1 
1 1 



.T 


n + 1 
1 2 




1 


. . . (B.3) 


The power output at x = 1 (hot junction), as derived in 

section 4.3.2 is given by : 


Ct) = .T^t) . ICt) + . 


aT^ (X , t) 
dx 




aT^cx, t) 


x= 1 


ax 


' x= 1 

... (4. 17) 


discretising it we have : 


Pi = - ^-^1 


/- n + l _ „n+l 

Mm Mm-1 
Ax 


+ k^.a^. 


r + 1 

Mm 2M-1 

Ax 


where T denotes the temperature of the first element at the 
1 M 

point ^iv6n by x = 1 , and siniilarly for th© other t errUd . 

From the boundary condition the temperatures of the two 
elements at x = 1 must be equal for all times, so that 
T = T 

IM 2M 

Hence the above expression can be rearranged as, 




Ax 


,T 


n + 1 
IM- 1 


^ M -^2 n +1 ^ 

2M- 1 
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.T 


n 4* 1 
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Designating the three matrices as A, T and B respectively the 
atxwe eq^)ation can be written in short form as : 

A.T = B 

The matrix T, ie. the temperature profile, is solved for on 
the conputer. Ar. exanple programme for the first case, le. the 
response of the thermoelanent for a step change in the current, is 
given in Appendix C. The prograimies for the other cases are very 
similar to the one listed in Appendix C. 
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APPENDIX - C 


PROGRAMME LISTING 




* 
* 
* 
* 
* 
* 
* 
* 

c 

Program Dynami c_Ana 1 ys i s_Of _A_Thermo_El ect ri c_Juncti on 

c 

c 

parameter (npara=200) 

Dimension x(npara), T(npara), v{npara ,npara) 

Real i , iO ,kl ,k2 , length 
common/ volt/n,m,dt, w_old 

common/ const/ al 0 , al 1 , yl ,y2 , 1 ength,ul ,u2 , al , a2 ,dx,plO ,p20 



c OPENING THE INPUT DATA FILE "in" 

open(ll,file='in' ) 



C OPENING THE OUTPUT DATA FILE 

open( 13 , f il e= ' stepO . 5 ' ) 



C READING THE INPUT PARAMETERS FROM THE FILE "in" 

read( 11 , * )p01 
read( 11 , * )dt , tmax 
read(ll,*)n 
read(ll , * )i0 
read( 11 , * )nitr 
m=2*n-2 




c THIS PROGRAMME CALCULATES THE HOT AND COLD JUNCTION 
c TEMPERATURES BY SOLVING THE GENERAL GOVERNING EQUATIONS 
c FOR THE TEMPERATURE DISTRIBUTION IN A THERMOELECTRIC 
c COOLING ELEMENT. THIS PROGRAMME IS FOR THE FIRST CASE 
c VIZ. A STEP CURRENT OF 0.5 A AND CONSTANT THERMAL POWER 
c INPUT OF 0.5 W. 
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ASSIGNING THE PROPERTY PARAMETERS OF THE THERMOELEMENT 

kl = 2.4 

k2=2 . 2 

length=2 . e-2 

ul= 0.91e-4 

u2= 0.91e-4 

wl=1.23e6 

w2=1.23e6 

yl = 0 .814e-5 

y2 = 0 .814e-5 

pl0=1.3e-5 

p20=l . 3e-5 

al = 0 .836e-4 

a2 = 0 . 836e-4 

al0=2.3e-4 

al 1=2 . 3e-4 

dx= 1 ength/m 

T0=25.0 

p0 = 0. 

pl=0. 
w_old=0 . 


INITIALIZING FOR TEMPERATURE 
do j=l , (m+1 ) 

T( j)=TO 
enddo 


time=0 
itr = 0 
i = 0 . 

write(13,*) 'Temp, profile for a step current of 0.5 amperes' 
write{13,*)' 

write(13,*)' TIME HOT JUNG. COLD JUNG. ' 

itr=itr+l 

time=time+dt 


1000 
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if (time.ge.5) go to 2 
go to 3 

2 i = iO 

p0=p01 

pl=-p01 + volt*i +w_rate 

3 continue 



c CALCULATING THE COEFFICIENTS OF THE DISCRETIZED 

c EQUATIONS 

c 

c COEFFICIENTS OF EQUATION B.l 

bl = ((kl/dx)+(ul*i)/(2*al))*(dt/(wl*dx)) 
b 2 = ( (-2*kl)/(wl*dx**2) + (yl*i** 2 )/(wl*al** 2 ))*dt + 1 . 
b3 = ( (kl/dx)-(ul*i)/(2*al))*dt/{wl*dx) 
b4 = ( (pl0*(i**2)*dt)/(wl*al**2)) 



c COEFFICIENTS OF EQUATION B.2 

dl = ((k2/dx)-(u2*i)/(2*a2))*dt/{w2*dx) 
d2 = ((-2*k2)/(w2*dx**2)+(y2*i**2)/(w2*a2**2))*dt +1. 
d3 = ((k2/dx)+(u2*i)/(2*a2))*dt/(w2*dx) 
d4 = ( (p20*(i**2)*dt)/(w2*a2**2) ) 


COEFFICIENTS OF EQUATION B.3 
cl = ((kl*al)/dx-(al0*i)+(k2*a2)/dx) 
c2 = -((kl*al)/dx) 
c3 = - ( ( k2*a2 ) /dx) 


COEFFICIENTS OF EQUATION B.4 
el = -((k2*a2)/dx) 

e2 = (((kl*al)/dx)+(all*i)+((k2*a2)/dx)) 
e3 = -((kl*al)/dx)' 


ASSIGNING THE TERMS TO THE R.H.S. MATRIX 
x( 1 )=p0 
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do j=2,n-l 

x( j)=bl*T( j + l)+b2*T( j)+b3*T( j-l)+b4 

enddo 




x(n) =pl 



do j=n+l,m 

x( j)=dl*T( j + l)+d2*T( j)+d3*T( j-l)+d4 

enddo 



c ASSIGNING THE TERMS TO THE COEFFICIENT MATRIX 

do j=l,m 
do k=l,in 
v( j,k)=0.0 
enddo 
enddo 



v(l,l)=cl 

v(l,2)=c2 

v(l,n+l)=c3 

v(n,in)=el 

v(n,n)=e2 

v(n,n-l)=e3 



do j=2,n-l 
v( j, j)=1.0 
enddo 




c 


do j=n+l,m 
v( j, j)=1.0 
enddo 


call gass(v, x,T,m) 

call vol tage( i , t , vol t , w_rate) 

if( mod(itr , nitr ) .eg. 0) then 
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wri te(13 , 222 ) time , T( 1 ) , T(n) 
print*, time, T(l) ,T(n) 
endif 

if( time .It. tmax) go to 1000 
1 continue 

rewind 13 

222 format(3fl5.5) 

stop 
end 



c THIS SUBROUTINE CALCULATES THE VOLTAGE BY SOLVING EQUATION * 

c 4.25, USING THE TEMPERATURE VALUES CALCULATED IN THE * 

c PROGRAMME . * 

c * 



subroutine VOLTAGE(i , t , vol t , w_rate) 
parameter (npara=200) 

Dimension T(npara) 
common/ VO 1 t/n ,m, dt ,w_old 

common/ const/alO ,al 1 , yl , y2 , 1 ength,ul , u2 , al , a2 ,dx , plO , p20 

real i , 1 ength 

suml=0 . 0 

sum2=0 . 0 

do k=l,n-l 

tav=0 . 5* ( t (k ) +t (k+1 ) ) 

suml=suml+tav*dx 

enddo 

do k=n+l,m 
sum2=sum2+ t(k)*dx 
enddo 

sum2=sum2+ 0 . 5*dx* ( t (n ) +t ( 1 ) ) 
w_new=al*suml + a2*sum2 
w_rate= ( w_new-w_ol d)/ dt 
w_old=w_new 

terml= al 1 *t (n) -al0*t { 1 ) 



102 


term2= -i*(length) * (pl0/al+p20/a2) 

term3= -i*yl*suml/al 

term4= i*y2*sum2/a2 

vol t=terml+term2+term3+term4 

return 

end 



c THIS SUBROUTINE SOLVES FOR THE TEMPERATURE DISTRIBUTION * 

c IN THE ELEMENT BY USING THE GAUSS-JORDAN ELIMINATION * 

c TECHNIQUE. * 


SUBROUTINE GASS ( A , Y , XX , N ) 
parameter (npara=200) 

DIMENSION A(npara ,npara ) , Y(npara) , XX(npara) 
M=N-1 

DO 10 1=1, M 
L=I + 1 

DO 10 J=L,N 
IF(A(J,I))6,10,6 
6 DO 8 K=L,N 

8 A(J,K)=A( J,K)-A(I,K)*A(J,I)/A(I,I) 

Y(J)=Y( J)-Y(I)*A(J,I)/A(I,I) 

10 CONTINUE 

\ XX(N)=Y(N)/A(N,N) 

DO 30 1=1, M 
K=N-I 
L = K+1 

DO 20 J=L,N 

20 Y(K)=Y(K)-XX( J)*A(K, J) 

XX(K)=Y(K)/A(K,K) 
iO CONTINUE 

RETURN 
END 



